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In this work, we calculate radiative vertex VPy (V = cj>,ui,p and P = 77,77') by utilizing ui — <j> 
mixing scheme and taking into account the contributions of the three-particle twist-4 distribution 
amplitudes of the photon in QCD sum rules on light cone. According to experimental data of 
V — > P7 and P — s> V7 from PDG, a value of the Tj — rj mixing angle, tp — (40.9 ± 0.5)°, is extracted 
in the framework of the quark-flavor basis to describe the 77 — 77' system. 
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I. INTRODUCTION 

The problem of 77 — if mixing in the pseudoscalar-meson nonet has been studied many times in the last forty years. 
Because those researches play a important role for the 5i7(3)-breaking effect and the U(1)a anomaly. In contrast to 
the 4> and uj mesons, where they are taken as almost ideally mixed states with quark content of well defined flavor, 
rj — if mixing is still a debated subject. In the pioneering work[l|-[3|, a mixing angle 9p was conventionally introduced 
to describe 77 and vf as linear combinations of octet and singlet basis states. With the development of the experiment, 
a phenomenological investigation [|| found that one mixing angel was insufficient to describe more physical processes, 
where 77 or rf meson appears in the initial state or the final state. And then two new equivalent schemes to describe 
the 77 — 77' mixing was proposed by LeutwylerQ and Feldmann ct al.Q, respectively. The correspondence mixing 
angels arc 9±, 9% for the octet-singlet basis and ^ q ,^p s for the quark-flavor basis. More literatures about two mixing 
angles can be found in Rcfs.[8Hl2|. where 9\ is not equal 9s, but <p q is equal ip s apart from terms which violate the 
Okubo-Zweig-Iizuka (OZI) rule[7[. In our next calculation, we chose ip q — cp s = Lp in the quark-flavor scheme to 
investigate radiative vertex VPj (V — 4>,to, p and P = 77, 77'). 

The radiative decays between light pseudoscalar (P) and light vector (V) mesons are an excellent laboratory for 
investigating the nature and extracting the non-perturbative parameters of light pseudoscalar nonet in low-energy 
hadron physics. Among the characteristics of the electromagnetic interaction processes, the coupling constant, <7yp 7 , 
plays one of the most important roles, since they determine the strength of the hadron interactions. So by investigating 
the above six radiative vertexes, in which 77 or 77' meson was involved, one can extracting the mixing angle of 77 — 77' 
system. Interest on this issue has been performed by many authors, for example, in Rcfs. (5, foL [T3T - [l6j . In this paper, 
we renewed this issue in QCD sum rules on light cone. Since the radiative decays of light meson is belong to the 
low-energy hadron interaction, which is governed by non-perturbative QCD, it is very difficult to obtain the numerical 
values of the coupling constants from the first principles. In order to interpret coupling constants from the experimental 
data, we immediately need to deal with large distance effects from the photon besides the hadrons, because a special 
feature of the QCD description of hard exclusive processes involving photon emission is that a real photon contains 
both a hard electromagnetic and a soft hadronic component. In R ef . njj , [l8j . a consistent technique was proposed by 
closed analogy with distribution amplitudes (DAs) of mesons [l9|, |2C|. The soft hadronic components of the photon 
are related to matrix elements of light-cone operators with different twist in the electromagnetic background field and 
can be parameterized in terms of photon DAs. Since the photon emission from the light quark takes place at large 
distances, the use of standard QCD sum rules based on the local operator product expansion (OPE)is not sufficient. 
Rather, one should use a light-cone expansion which is adequate for exclusive processes with light particles. 

The paper is organized as follows. In section 2, we derive the coupling constants gvp-y by utilizing uj — <p mixing 
scheme and taking into account the contributions of the three-particle twist-4 distribution amplitudes of the photon in 
the light-cone sum rules. In section 3, wc present our numerical analysis. The final section is reserved for a conclusion. 
The photon distribution amplitudes arc list in Appendix A and the overlap amplitudes for pseudoscalar mesons and 
vector mesons, which were defined in section 2, are presented in Appendix B. 
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II. RADIATIVE VERTEX VPj IN LIGHT-CONE SUM RULES 

In the framework of light cone QCD sum rules, we immediately choose the two point correlation function with the 
photon as follow 

T,{ Pl q) = i [ d^xe-^ < 7 (p)|T{,r(x)j 5 p (0)}|0 >, (1) 



to extract the radiative coupling constant gvp-y- Here jY(x) is the vector meson current and jrf (0) is the pseudoscalar 
current. According to the quark model of hadron and neglecting the contribution from the high Fock state of hadron, 
the above interpolating currents may be written as 

jf = + si 11 / 3 + sj^s cos /3, j£ = -^(uj^u + d-f^d) cos /3 - s-f^s sin (3, j? = -^=(u-f^u - d-f^d) (2) 

for light vector mesons 0, w, p and 

jg = — y=(uijsU + di'ysd) cos (p — si'jsssmip, = —=(uij5U + dij$d) simp + sij^s cos/p (3) 
v 2 v 2 

for light pseudoscalar mesons 77,77' , respectively. Here <p is the value of the 77 — 77' mixing angle, which will be discussed 
in this work, and (3 is the value of the w — <fi mixing angle, which has been determined from the available experimental 
data in the Ref.[2l| as /3 = 3.18°. 

According to the basic assumption of quark-hadron duality in the QCD sum rules approach, we can insert two 
complete series of intermediate states with the same quantum numbers into the correlation function to obtain the 
hadronic representation. After isolating the contribution of the ground state by the pole terms of the vector meson 
and the pseudoscalar meson, we get the following result 

rp , v _ < P(-qh(p)\v(p- q) >< v(p - g)Lff |Q >< QbTl^-g) > ... 

W'^" {{p _ q) 2_ m 2 v){{ _ q? _ m 2 p) +•••> W 

for V —> P7 decay, where —q, p and p — q denote the pseudoscalar meson, the photon and the vector meson momentum, 
respectively. The amplitudes of these interpolating currents with the meson states are defined as 

<0\jV\V(p-q)> = X v ul (5) 
< 0\£\P(-q) > =Ap, (6) 

where is the polarization vector of the vector meson, Ay and Ap are called the overlap amplitudes which can be 
determined by QCD sum rules method in Appendix B. The coupling constant <?yp 7 is defined through the effective 
Lagrangian 

2 = -—gvp^^id^^y - d v ^)(d a A fi - dPA a )<t> P (7) 
my 

where 0y,0p and A denote the vector field, the pseudoscalar field and the photon field, respectively. Therefore, the 
< P(— q)^(p)\V(p — q) > matrix in the hadronic representation can be written as 

< P(-qh(p)\V( P - q) >= -^-g V P 7 K(p 2 )s^ulq^lp^ (8) 



where p 2 — for the momentum of the real photon and K{p 2 ) is a form factor with K(0) = 1. Substituting eq.(5),(6) 

-vafi _7 



and (8) into the hadronic representation, we obtain the physical part and choose the structure e u ""p y q a el from which 



the corresponding invariant amplitude, 

T{{p - q) 2 , q 2 ) = j- 6 ^2 PlX 2U 2 2T + / 7 9 {yS }^ 2)ds ^ dS2 + subtractions. (9) 

mv {(P - q) z - ) {q z - mf,) 7 S (si - q l ) (s 2 - (p - q) 2 ) 

The first term is the contribution of the ground-state and contains the gyp-y coupling, while the hadronic spectral 
function p h (s\,S2) represents the contribution of higher resonances and continuum states. The integration region in 
the (si, S2) - plane is denoted by £ and one may take (si) a + (s2) a < ( s o) a , where sq is the effective threshold in the 
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double dispersion relation. It is relevant with Sq and s^ , which are the effective thresholds in the vector meson and 
the pseudoscalar meson channels, respectively. How to take their values will be discussed in the next section. The 
subtraction terms isn't considered in the sum rules, because they will be removed by a double Borel transformation. 

Next, we calculate the correlation function from QCD side by using light cone operator product expansion method, 
in which we work with large momenta, i.e., —q 2 and — (p — q) 2 are both large. The correlation function, then, can 
be calculated as an expansion near to the light cone x 2 w 0. The expansion involves matrix elements of the nonlocal 
operators between vacuum and the photon states in terms of the photon DAs with increasing twist. At the same 
time, the full quark propagator of the light quark (22l . |23| in the presents of gluonic and electromagnetic background 
fields is used in this calculation, and it is given as 



iS{x,0) = < 0\T{q(x)q{0)}\0 >= i 



2ir 2 x 4 4ir 2 x 2 



<qq> (l 1 i mq ^ 



12 



m n < qq > (1 + i 

192 m 1 6 



1 i 

ig s j dv{ 
o 

i 

dv{ 



16n 2 x 2 



Til 

iyto^ + va^)G^{vx) ~ i—^G^{vx)a^ ln(- 



o 



2 7 s)} 



H — - A — + 2 1e)} 



(10) 



where G M „ is the gluon field strength tensor and F^ u is the electromagnetic field strength tensor. Taking into account 
the flavor SU(3)-breaking effect, we take the quark mass as m u = rru = and m s ^ in our calculation. 

We substitute the above propagator and the relevant photon DAs Qjillll into correlation function (1) and integrate 
over the time-space coordinate x, we can obtain the expression of the amplitude T((p— q) 2 , q 2 ) up to twist-4 accuracy. 
The following taking place is a standard process in light-cone sum rules. The hadronic spectral function p (si, S2), 
which was used to control the contributions of excited states and of the continuum, can be approximately evaluated by 
using quark-hadron duality. To suppress the contributions of the excited and continuum states and of the subtraction 
terms in the hadronic representation, the double Borel transformation [24T - |26j about the variables — q 2 and — (p — q) 2 
was introduced and the useful formulas for transformation as follow: 



T(a) 



uq 2 — u(j> — q) 2 



(M 2 ) 2 - a u 
-MfMl 5{u - Uo) 



(11) 



and 
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Here M 2 = jjr—jfp-, uq = M -i+ M -i , M'f and Mf arc the Borel parameters associated with — q 2 and — (p — q) 2 , 
respectively. 

After the above lengthy calculation, we obtained the final result for the coupling gvPj- 



1 {{p-q) 2 - m 2 v ){q 2 - m 2 p )\ M 2 M 2 



1 -^p. 

e M i e M £ . 



(12) 



1 ^v+^-P 
AyAp 



A VP <qq> 



X + B VP <ss> (X + B) 



(13) 



where X = A + Ig(uq) + If(uq), 



-M Z X 



k 



/feV (l-2ttp) 
y 2> Z-« r (k-n)\n\ 

n=0 v ' 



k- n k—n 1 so \j 



A(uq) 



(14) 



B = 



_- M ^(l-e-^)(l-^) 



-^W°>(uo) 



2 < ss > 2 < ss > 



ZM 2 



-mi + 



< ss > 
3tt 2 M 4 ' 



(15) 



IgM 



Uo rl—uo 

da„ / da„- 



1 







1 1 

1 — Ctq - Uq 

U pl—UQ 

S(an))]-2 I da q 

Jo 



[h{o.i) - T 2 {a>i) + T 3 (ai) - 7i («») - <$(«») 



da <i7T — — — -\^[Ts(oh) - %(ai) - S(on)) 
(1 - a q - a q y 



(16) 
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and 

I f (uq) = - / da q / dctq- [5 7 (ai) + 77(aO] + 2 / da q \ da q — T?77(ai)- ( 17 ) 

Jo Jo 1- a q - a q J J (1 - a q - a q y 

Here x is t ne magnetic susceptibility of the quark condensate, which has been introduced in the pioneering work [27j for 
proton and neutron magnetic moments. In our next numerical analysis, we take xi^GeV 2 ) = — 3.15±0.3Gel/ -2 which 
was obtained by using QCD sum rules analysis of two-points correlation function [I?)- fj{u) is the photon leading 
twist distribution amplitude, 0^°' (u ) is the photon twist-3 distribution amplitude, A(u), Tj(ai)(j — 1, 2, 3, 4), 7^(ai), 
5(ai), S(oti), and S 7 (ai)(i = q,q,g and a g = 1 — a q — a g ) are the photon twist-4 DAs. Their detailed expression 
are given in Appendix A. = — (0.0039 ± 0.0020)GeI/ 2 [l7( is the nonperturbative constant to describe the photon 
twist-3 DAs and 7e is Eulcr constant. bj~ are the coefficients of the leading twist-2 distribution amplitude (p-y(u) by 
exploited as a power series in (1 — u), (p-y(u) = J2k &fc(l ~ u ) k ■ The coupling gvp-y for the P —> Vj decay can also be 
calculated by the similar approach. A summary of the results is presented in Table 1 . 



\. p 


V 


v' 


A V p 


Bvp 


A V p 


Bvp 


<t> 


(e u + e d ) sin j3 cos ip 


—m^e s cos P sin ip 


m ( f 1 (e u + e d ) sin j3 sin tp 


m^es cos ft cos tp 




(e„ + e^) cos /3 cos <p> 


m„e s sin f3 sin ip 


m v > (e„ + e d ) cos /3 sin ip 


—m^i e s sin j3 cos ip 


P 


m p (e u - e d ) cos ip 





m n i(e u — e d ) smtp 






Table l:The parameters for the coupling gypj in the radiative decays of the light mesons. 



III. NUMERICAL ANALYSIS 



Now, we present our numerical analysis of the coupling constants gypj for the V —> Pj and P — s- Vj decays. 
In order to obtain numerical results of the sum rules from Eq.(13) and Table. 1, we take the input parameters as 
usual: m = 1.02GeV, m w = 0.782GeF, m p = 0.77GeV, m, ; = 0.55G e y, Try = 0.958GeF, m s = 0.156GeF, (qq) = 
-(0.24GeV r ) 3 , (ss) = 0.8(qq). The values of overlap amplitudes Ay from Eq.(5) are A = 0.250 ± 0.009GeV 2 ,A w = 
0.162 ± 0.00406^, A p = 0.150±0.003GeV^ 2 and their numerical analysis was presented in Appendix B. There haven't 
the numerical results for the overlap amplitudes Ap from Eq.(6), but their expressions with the mixing angle ip were 
given by using QCD sum rules in Appendix B. 

The following issue is to discuss how to reasonably choose the effective thresholds so,Sq,Sq for the sum rules 
of the coupling gvp-y- Sq and Sq are the effective thresholds in the vector meson and the pseudoscalar meson 
channels, respectively. In general, their values ranges from the mass square of the ground state to the mass square 
of the first exciting state. According to data of PDG HH , we can find the first exciting states of the above meson. 
They are </>(1680) for 4> meson, w(1420) for lj meson, p(1450) for p meson, ?y(1295) for r\ meson, r/(1405) for rj' 
meson. While taking into account the especial property of the double spectral function at si = s 2 , we obtain the 
eventual range of the effective threshold in the double dispersion relation: sq <G (1-04, 1.68)Gey 2 for the coupling 
g^-n-y, s G (1.04, 1.97)Ge^ 2 for the coupling g^, s G (0.61, \M)GeV 2 for the coupling g um , s G (0.59, 1.68)Gel^ 2 
for the coupling g pr)1 ^ sq G (0.92, 1.97)GeV^ 2 for the coupling g uv '~f and the coupling g pn '-y. In the next analysis, the 
effective threshold so will be strictly choose in the above ranges to fit the result of the sum rules and experimental 
data. 

In Fig.l, we discussed the Borel window of our sum rules. Here the mixing angle (p is fixed at ip = 40°, the Borcl 
parameter Mf of every coupling constant takes a defined value and the effective thresholds sq of the decay channel 
have different value which belong to the above discussing ranges. We find that there have a platform, where the 
coupling gvp-y is practically independent of the Borcl parameter M 2 , 1.8GeV^ 2 < M 2 < 2.8GeV 2 for the r\ channel 
and X.hGeV 2 < M 2 < 2.5GeV 2 for the rf channel. So our sum rules are reasonable and significative. 

In Fig. 2, we show theoretical and experimental values of the couplings gyp-y as functions of the r\ — ?y' mixing angle 
ip with the defined so, M 2 and the variable M 2 . The ranges of M 2 are the Borcl windows from Fig.l. So theoretical 
predictions are presented in the shadows in Fig. 2. By using the definition of gyp-y from Eq.(8), the decay widths of 
V — > Pj and P — > V~f are written as 

r(F^p 7 ) = a t P7(< " 5 m2p)3 . nP^Vi) = a9 \ pAm2p -f v) \ (18) 

24 my 8 nip 
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where a = 1/137 is the electromagnetic coupling constant. Comparing Eq.(18) with experimental data [28[, we obtain 
the values of the coupling constants, which are presented in the dot-dashed curves in Fig. 2. Finally, the pseudoscalar 
mixing angle <p from different V — rj, rf electromagnetic coupling processes and their average value are list in Table. 2. 
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FIG. 1: The the coupling constants gvp-y and gpvj as a function of the Borel parameter M ± for different values of the 
threshold parameters s with defined Af|. (a) M| = 2.0GeV 2 ; (b) Mf = 2.3GeV 2 ; (c) Mf = 1.9Ge1/ 2 ; (d) M| = 2.1GeV 2 ; 
(e) M| = 1.7GeV 2 ; (f) Mf = 1.9GeV 2 . 
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FIG. 2: Theoretical and experimental values of the couplings gvp-y as functions of the rj — rj' mixing angle ip. The shadows 
are theoretical predictions according to Eq.(13) and the dot-dashed curves are the couplings extracted from experimental data 
[H by Eq.(18). (a) Mf = 2.0GeV 2 and s = 1.5GeV 2 ; (b) Mf = 2.3GeV 2 and s = l.6GeV 2 ; (c)M| = 1.9GeV 2 and 
s = 1.5Ge1/ 2 ; (d) Mf = 2.1GeV 2 and s = 1.5GeV 2 ; (e) Mf = 1.7GeV 2 and s = 1.5GeV 2 ; and (f) M 2 2 = 1.9GeV 2 and 
so = 1.2GeV 2 . 



V(P) -> P(V) 7 


— > 777 


w — > 777 


p -»■ 777 


— ^ 77'7 


77' — ► 017 


77' -> p7 


<y3au(°) 


v(°) 


40.2 ±0.7 


41.2 ± 5.4 


40.3 ±3.4 


40.8 ± 1.2 


41.6 ±2.4 


41.1 ± 0.4 


40.9 ±0.5 



Table 2:The mixing angle ip from different V — 77, 77' electromagnetic coupling processes. 
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IV. CONCLUSIONS 



We calculated the coupling <?yp 7 (V = 4>,uj 7 p and P = r/ 7 if) of the V — > P7 and P — > Vj electromagnetic decays in 
the light-cone QCD sum rules. Comparing theoretical results and experimental data, we extracted a new pseudoscalar 
mixing angle p> = 40.9 ± 0.5° in the quark-flavor basis. This result is in agreement with Ref.Q, where the average 
if = 39.3 ± 1.0°. Recently, the KLOE Collaboration [H HI has measured the ratio i? = B{(p -> rfi)/B(^> ->• r/7) = 
4.77 x 10~ 3 , the pseudoscalar mixing angle tp = 41.4 ± 1.6° with the zero gluonium content for 77' and (p = 39. 7± 0.7° 
with the gluonium content for 77'. There is a little discrepancy between our theoretical results, R$ = 4.85 x 10 -3 and 
tp = 40.5 ± 1.0° from the <p — > 7/7 and <f> — > 7/7 decays, and the experimental results from KLOE. A possible reason is 
that we should consider the contribution from the gluonium content of 77' meson in our calculation. This will be our 
next work. 



Appendix A: The photon distribution amplitudes 

In this section, the clear expressions for the photon distribution amplitudes are showed as[l7l Il8j 

tp-y{u) = 6uu{l + 1P2C2 (u — u)) , 

^\u) = (l-(27 i -l) 2 )(5(27 i -l) 2 -l)^(l + A<-l^), 

A(u) = 40uV(3k + + 8(^-3(2) [ww(2 + 13u7t) 

+ 2it 3 (10 - 15u + 6u 2 )ln(u) + 2-u 3 (10 - 15t* + 6u 2 )ln(u)] , 
T\(ati) = -120(3(2 + Ct)( a g - a q )a q a q a g , 

T 2 (a 4 ) = 30a2(a 9 --a 9 )(( K - K + ) + (Ci-Ci + )(l-2a s ) + C2(3-4a 5 )), 
75 (a») = -120(3C 2 - C2" ~ Oi q )a q a q a g , 

T 4 (a 4 ) = 30a2(a 9 --a 9 )(( K + K + ) + (Ci+Ci + )(l-2a g ) + C 2 (3-4a ff )), 
5(«i) = 30a 2 g {(K + K+)(l-a g ) + (( 1 +(+){l-a g )(l-2a g ) 

+ C2[3(a 9 - - a q ) 2 - %(1 - a g )]}, 
S{ ai ) = -30a2{( K - K +)(l- ag ) + (C 1 -C 1 + )(l-a g )(l-2a 9 ) 
+ C2[3(a 9 - - a q ) 2 - a g (l - a g )]}, 
Sj(ai) = 60a 2 (a q + a,)(4 - 7(a q + a,)), 
77(«i) = 60ag(a 9 - a§)(4 - 7(a q + a q )). 

Here tpj(u) is the photon leading twist distribution amplitude, (uq) is the photon twist-3 DA, A(u), 7j(cti)(j = 
1,2,3,4), T^ipti), S{a,i), S(ai), and <S 7 (q:,-)(« = q,q,g and a ff = 1 — a q — aig) are the photon twist 4 DAs. The 
parameters appearing in the above DAs are given as pi = 0, n = 0.2, k + = 0, Ci = 0.4, uY = 3.8, uj^ = —2.1, 
C2 = 0.3, Ci + = and ^ = at the scale p = lGeV. 



Appendix B: Calculation of Overlap Amplitude 



In this section, we present the discussion of the overlap amplitudes for the vector mesons and pseudoscalar mesons. 
They are important input parameters for the coupling gyp-y in our sum rules. The overlap amplitudes A r) , A^' for 



pseudoscalar meson 77 and 77' are given as 31 



A, 



A^ cos ip — \ s n sin tp, X v i = A*, cos p + X^, sin tp 



in QCD sum rules, where 

(X q P ) 2 = e^M 4 { 



(Ap) 2 



e m- 



M 4 { 



87T 2 A/ 2 



[1-(1 



8M 4 
8AI 4 



1127T 



M 4 27M 6 



m s (ss)) 



1 



1127T 



M 4 27M e 



< ^/a~ s qq > 2 }, 

< sfalss > 2 }. 
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Here Sq is the effective threshoid and the pseudoscalar mixing angie <p will be determined by numerical analysis of 

With similar processes mentioned in the overlap amplitudes Xp of pseudoscalar meson, the overlap amplitudes Ay 
of vector mesons can be obtained and Xy S are written as prjj : 



(A y ) 2 



(Ay) 2 



mlM 2 e m -' M2 



< ^G%G a ^ > H27T a s < uu > 5 



12M 4 

m 2 v M 2 e m v' M " 



81 



36A/ 8 



A 



4^ (1 " e ° 7 >< 1 + T )+ A/* <M> 



< ^-G% v G a ^ v > 112tt a s < ss > 2 



12M 4 



81 



M 6 



18M 8 



<9sSa^^G a ^s> 



And then we take the input parameters: — 1.02GeV, m w = 0.782GeV, to p = 0.77GeV, m u = 0.005GeV, md = 
0.008GeV", m s = 0.156GeV, (qq) = -(0.24GeV) 3 , (ss) = 0.8(qq), a s = 0.5, m 2 , = 0.8GeV 2 , ^G 2 = 0.012GeV" 4 , 
< 9sS<7 li yk^-G a ^ 1 ' s >= m\ < ss > into the above expressions and yield: 

A = X q , sin/3 + X s , cos (3 = 0.250 ± 0.009GeV" 2 , 



A w = A« cos/3 - A^ sin/3 = 0.162 ± 0.004G e y 2 



X p = X q p = 0.150 ± 0.003GeV 2 . 
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